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Abstract 

We calculated numerically the localization length of one-dimensional Anderson 
model with diagonal disorder. For weak disorder, we showed that the localiza- 
tion length changes continuously as the energy changes from the band center 
to the boundary of the anomalous region near the band edge. We found that 
all the localization lengths for different disorder strengths and different energies 
collapse onto a single curve, which can be fitted by a simple equation. Thus 
the description of the perturbation theory and the band center anomaly were 
unified into this equation. 
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1. Introduction 

Electronic transport properties, the motion of electrons, in a random poten- 
tial are closely related to the phenomenon of localization. Since the pioneering 
work of Anderson Q for disordered systems fifty years ago, Anderson localiza- 



tion has been applied to various fields including photonics and cold atoms [15| . 
However, even for the simplest one-dimensional case, we have not found exact 
analytical treatment for arbitrary disorder strengths and energies 14 , Q ■ 



Accurate numerical approaches have been developed by using the quantum 
transfer matrix rcnormalization group method for finite temperature systems 
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25 1, the density matrix renormalization group method for interacting systems 
20 1 , and the integral equation method for systems in the thermodynamic limit 
11], respectively. In this work we will give a unified description of the in-band 
(without the band edge anomalous region) behavior for the one-dimensional 
model with weak diagonal disorder. We study the problem with the help of a 



numerical method we developed earlier in Ref . , which was an application of 
the transfer matrix method jl8l | in localized phase in the thermodynamic limit. 

The localization theory shows rigorously that all the eigenstates are expo- 
nentially localized for one-dimensional uncorrelated disordered systems. The 
single parameter scaling (SPS) theory [l[ gave great insight into the properties 
of disordered systems. It argues that the dimensionless conductance g is the 
only relevant parameter that controls its variation with system size L. SPS re- 
quires that the complete distribution function p(g) should be determined by one 
single parameter. Originally, SPS was derived within the random phase approx- 
imation [3[. The finite size Lyapunov exponent j(L) = (l/2L)ln(l + l/g) was 
proposed to be an appropriate scaling parameter (it approaches the nonrandom 
limit when L — > oo [17j]). It is well-known that the random phase approximation 
fails at the band edges E = ±2 as well as the band center E = 0, H3|. An- 
other length scale was suggested to give an appropriate criterion for the failure 
of SPS [p and a relatively good scaling parameter was found for the non-SPS 
region [6|. 

For weak disorder, although p(g) has two independent parameters which is 
in contrast with the case of strong disorder, it was found that one of them is a 
universal number [lij . Thus it is still SPS. There exists a perturbation theory 
in which the Lyapunov exponent depends on energy E and disorder strength W 



24] 



^ W2 )= 96(1^/4) ' W 
At the band center, a revised perturbation gives [l2|, [H, E3] 

W 2 

1 (W 2 ) = . (2) 

n ' 105.045- •• v ; 

There was no analytical expression which connects smoothly the above two 
perturbative results. To clarify the smooth crossover in between we study the 
Lyapunov exponent with two variables, the energy E and the disorder strength 
cr in this paper (a 2 — W 2 /12). Furthermore, we only study the problem in the 
thermodynamic limit and thus the third variable, the size of a finite system, is 
not in consideration. 

It was shown that for one-dimensional Anderson model with diagonal disor- 
der, the band center is actually a band boundary adjoining two neighbor bands 
rather than the center of a single band, which is the reason for the fail of SPS 
near the band center. Therefore it is similar to the band edge 0. This should 
also be why in Ref. [10| the authors used a four-step map at the band center. 
In our recent work, we proposed a parametrization method of the transfer ma- 
trix for one-dimensional Anderson model with diagonal disorder. Making use of 
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this method, we obtained numerical results of the dependence of the Lyapunov 
exponent on energy and disorder strength in the thermodynamic limit within 
the localization regime. A remarkable coincidence of our numerical results with 
the existed results demonstrated the reliability and efficiency of our method 



For the anomalous region near the band center, there have been analytical 



results for E — >• and E ^> a 2 through the scaled energy E/a 2 [5|. In the 
present work, we will show the detailed numerical results about the anomalous 
behavior of the Lyapunov exponent with weak disorder near the band center 
by using our method, which coincides with the known analytical results; and 
we will give a simple equation to describe the behavior for arbitrary weak dis- 
order strength and for all E in the band except the anomalous region near the 
band edge through the parameter E/a 2 . Thus the band center anomaly and 
the perturbation theory are connected smoothly. 



2. Parametrization method 

Consider the one-dimensional Anderson model with diagonal disorder, 

ipi-x + ip i+ l = (E - ei)ipi, (3) 

where ipi is the electron wavefunction at site-i and is the on-site energy from 
a random distribution p € {e). For the box distribution, p € {e) = 1/W, \e\ < W/2, 
(e 2 ) = W 2 1 12, and for the Gaussian distribution, p £ (e) = (l/v^cr) exp(-e 2 /2<r 2 ), 
(e 2 ) = a 2 . In the transfer matrix method, Eq. ^ can be written as 

where Vi = E — ej and Tj is the transfer matrix. 

Using a parametrization method of the transfer matrix proposed in our pre- 
vious work [lll | , we can calculate the Lyapunov exponent in the thermodynamic 
limit within the localization regime. Let = TlT^_i • • • TV Then we can 
parameterize MM* as follows 

U(^)M x MiU(-0 L ) = ( ^ £ _ Xl \ , (5) 

where M* is the transpose of M and 

U(0 L ) = ( C °1 L " Sm / L V (6) 
v ' V sm6>£ cos6>l / 



The recursion relation of 6 in large L limit is 



tan0 i+ i = 1 . (7) 

vl+1 - tan&L 
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The two parameters A and 9 are what we need in order to calculate the Lyapunov 
exponent 7. The equations we obtained for the distribution function p(9) and 
the Lyapunov exponent are 

= -r?Z I P(0')pv(^+ttm6')d9' (8) 



sm 



7 = - J p(6)\n\tan8\d8 (9) 

where p v (v) is the random distribution function of v. In this study we use the 
Gaussian distribution p v (v) = (l/V^nr) cxp[— (v — E) 2 /2a 2 ] to solve Eq. ® 
and to calculate 7 numerically for different a 2 and different energy E. 



3. Numerical results 

Fig. [Hand Fig. [2] illustrate the distribution functions p(6) which were calcu- 
lated as p v (v) being the Gaussian distribution. In Fig. [U four curves of E = 
are plotted for different values of a 2 decreased from 5.3792 to 5.24288 x 10~ D ; 
and in Fig. [2l four curves of a 2 = 0.00125 are plotted for different values of E 
decreased from 0.8 to 0.05. In Fig. [lj p(0) becomes a stationary distribution 
when disorder is weak enough (the line of a 2 = 5.24288 x 10 -5 ) 

where K is the complete elliptic integral of the first kind [10(. There is a 
correspondence between p{9) and the distribution of the reflection phase <p in 
Ref. [4|, which gave similar plots, through the relation tan# = cos(4>/2)/ cos(/i + 
0/2), where 2 cos fi = E. In Fig. [21 the expression of distribution p(9) for the 
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Figure 2: (Color online) Distribution p(0) from Eq. JHJl at fixed <r 2 = 0.00125 for the Gaussian 
distribution with decreasing E. 



three larger values of E is 

= sin/z 

v ; tt(1 - cos/zsin2(9) v ' 

corresponding to the uniform distribution for the perturbation theory (l2l . [Toll . 
A recent numerical calculation [l3[ of the angle distribution function in Ref. [lfj| 
also gave a similar behavior. In our numerical approach we have set a relative 
precision 10~ 10 of p{9). 

For weak disorder we have known that the perturbation theory breaks down 
near the band center. It is clear in Fig. [3] that the distribution p(9) of E = 
0.05 with a 2 = 0.00125 deviates from Eq. ([TTj) . which means that the band 
center anomalous region extends up to E > 80a 2 . In Figs. [1] and [2] we have 
demonstrated known results for the two limits E/a 2 —> and E/a 2 — > 00; 
and in Fig. [3] contribution of higher order terms to the curves in weak disorder 
perturbation becomes apparent for finite E/a 2 . In Ref. [22| the authors gave 
a complete accurate analysis of not only the localization length, but all higher 
moments of the distribution of the Lyapunov exponent for finite systems. It 
is natural to expect that the deviation from Eq. (fTTj) comes from higher order 
terms in the perturbation given in Ref. [22I ] . And the present work will give an 
accurate equation for the localization length in the whole band except the band 
edge anomalous region. 

To demonstrate the anomalous behavior near the band center, we plot 'y(E) 
with disorder strength a 2 = 0.005 for p v (v) being the Gaussian distribution in 
Fig. 01 It shows clearly that in the vicinity of the band center the behavior 
of the Lyapunov exponent deviates from the perturbation result ("f p solid line) 
and when E approaches zero it becomes the result of Eq. @ smoothly [22j ]. 
Our numerical results coincides with the analytical results for the region near 
E/a 2 < 1 and E/a 2 > 1 [5]. We will present a quantitative description by 
collapsing the curves for different energies and disorder strengths through the 
parameter E/a 2 . 




Figure 3: (Color online) Comparison of p(8) between the numerical result and Eq. mil l (E = 
0.05 and a 2 = 0.00125), which means that the standard perturbation theory fails when E 
approaches <r. 
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Considering Eqs. JT]) and © as well as the result shown in Fig. |H we choose 
y = 1 — -f/'fp as a, scaling parameter, y is the relative difference between the 
perturbation result and the numerical result. To show all the values of E/cr 2 
we use x = arctan(£ , /acr 2 ), where a is a fitting coefficient, x = at the band 
center E — and x = tt/2 for a finite energy when E/a 2 ~ 1/c 2 — > oo with 
a 2 —> 0. In Fig. [5] there are three groups of data and the fitting coefficient 
a w 1.436. Each group has a fixed weak disorder strength. All the data points 
merge into a single curve. Thus the in-band behavior can be described by a 
single scaling parameter for all weak disorder values of a 2 as long as E is not 
located in the band edge anomalous region. If the energy is in the band edge 
anomalous region, we found that y deviates from the universal curve in Fig. [5j 
which was not shown here. 

When E approaches the boundary of the band center anomalous region, 7 
becomes equal to the perturbation result j p , therefore y — > 0; when E = 0, 
y = 1 - 96/105.045 w 0.086. So We expect that 

96 

^ (1 -ro5^45 )cos2 ^ (12) 

where cos 2 a; = [1 + {E/aa 2 ) 2 ]^ 1 . Therefore the Lyapunov exponent for weak 
disorder and with energy in the band except near the band edge is 

7 = /(^)7p, (13) 

where f(E/<r 2 ) = 1 - 0.086[1 + (E/acr 2 ) 2 }- 1 H . The difference in 7 is only 

8.6% for the E = and a finite energy. A quantitative description can be given: 
if we set the criterion for the breaking down of the perturbative result in Eq. (TTJ) 
as 1% difference in 7, then by Eq. (| X3[) we obtain \E\ > 4cr 2 as the region for 
Eq. (fTJ to be considered as valid. It is \E\ > 8D which is in agreement with the 
\E\ > 10D given in Ref. 22]. Therefore there is no definite boundary between 
the perturbation theory region and the band center anomalous region. 

In Fig. [5] all points collapse onto a single curve. We should mention that 
Eq. (|12[) only includes the main correction from E/a 2 . There are higher order 
corrections in f{E/a 2 ) The small deviation between Eq. (jT2"j) and the 

numerical data was shown in Fig. [51 

At x = point, where E = 0, it is already known that Eq. (|2|) is only valid in 
weak disorder limit. 7 behaves in logarithmic of the disorder strength in strong 
disorder limit, and the crossover curve between the two limits can be found in 



Ref. (11]. The difference around x = point is reasonable: a 2 = 0.02 is not 
weak enough and higher order corrections of the perturbation in a 2 becomes 
apparent. In the weak disorder limit {a 2 = 0.005, 0.00125 in Fig. [5]) all the 
curves for different a 2 is a single curve, which depends on E/a 2 only. 

4. Conclusion 

In summary, we calculated the inverse localization length in one-dimensional 
Anderson model with diagonal disorder. We obtained numerically the curve of 
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Figure 5: (Color online) Three groups of data fall on a single line with our choice of y and x. 
Eq,l ll2l) coincides well with the numerical data in this figure 
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the inverse localization length for energies inside the energy band in the case of 
weak disorder. A unifying curve was given for different weak disorder strengths 
and energies outside the band edge anomalous region. 

The scaled energy used to plot the unifying curve is E/a 2 , which has cor- 
respondence to the scaling parameter k = li oc /l s proposed by Deych et al. Q- 
They used k for the criterion of SPS, which has a clear physical image. However, 
the new length scale l s and the localization length li oc are both functions of the 
energy E and the disorder strength a 2 . Thus k is a two parameter function 
k(E,<j 2 ). And we showed in this work that k actually depends on the single 
parameter E /a 2 , which is the essential parameter to give the universal curve in 
Fig. 13 

We also found a similar unifying description of the perturbation theory and 
the band center anomaly for the box distribution of disorder. We believe that 
other kinds of uncorrelated disorder with a finite variance can have the same 
unifying description. It should be mentioned that the Lorentzian distribution 
has no such band center anomaly, because its variance does not exist. 
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